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Abstract 

^J^, We start with a formal series expansion of the regular asymptotic solution of the higher order 

gravity theory derived from the quadratic lagrangian R + eR^ in vacuum, and build a time 
^ development as a solution of an initial value formulation of the theory. We show that the 

^SJ ■ resulting solution, if analytic, contains the right number of free functions to qualify as a general 



solution of the theory. We further show that any regular analytic solution which satisfies the 
(N . constraints and evolution equations can be given in the form of such an asymptotic formal series 



(N 



expansion. 



1 Introduction 



X 

^ ■ Ever since the fundamental realization that Einstein's equations can be equivalently studied as a 

geometric, nonlinear system of evolution equations with constraints, there has been a considerable 
stream of problems in general relativity of a classical dynamical nature crucially depending on 
the notion of evolution from prescribed initial data (cf. [1] and references therein). In particular, 
for the so-called cosmological case, at least since the investigations of Lifshitz and Khalatnikov 
in relativistic cosmology [2], perturbative approaches to the question of genericity of cosmological 
solutions constructed asymptotically from given initial data, have been advanced in various direc- 
tions in general cosmological theory. This was the way to prove the existence of vacuum, regular 
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cosmological solutions in general relativity having the required number to qualify as general ones, 
and also the non existence of a general singular, radiation solution in the same context, cf. [21 E], 
eventually leading to the realization that the initial state was of a more complex nature [U [5l [6] . 
These first results were based on an approximation technique that consisted of several steps such as 
writing down a suitable expansion of the metric, substituting it to the field equations and counting 
the number of free (or arbitrary) functions needed to make the whole scheme consistent. 

The original perturbative approximation scheme mentioned above proved to be especially fertile, 
and in fact it is being exploited in various directions ever since. It was used by Starobinski to study 
solutions of the Einstein equations with a positive cosmological constant, especially with regard 
to the question of the asymptotic stability of de Sitter space used in inflationary scenarios [TJ [8] . 
The formal series expansions were used more recently in [9j to study the genericity question in 
relativistic cosmologies with a general equation of state p = wp, where it was proved that in 
order to be able to construct a general singular solution initially, certain restrictions on the fluid 
parameter w are needed. The ultrastiff case was recently considered in a more rigorous way using 
Fuchsian techniques in Ref. [lOj. Formal series where also used in [llj to study the more involved 
problem of perturbing an FRW universe containing two suitable fluids, and in [12] where it was 
shown how to construct a sudden singularity with the genericity properties of a general solution, a 
situation met previously only in the 'no-hair' behaviour of inflation. Extensions to certain higher 
order gravity theories have also been considered in a formal series context, especially in connection 
with the stability of de Sitter space under generic perturbations in these theories, cf. [13^ I14j. 

It is interesting that the original approximation scheme of formal perturbation series has been 
further refined and applied in various situations in cosmology, cf. [15^ \TE[ [17] and references 
therein. More recently, Rendall has put the original formal series expansion techniques used by 
Starobinski [7j for the stability of de Sitter space, in a more rigorous basis and was able to prove 
various interesting theorems concerning function counting and formal series solutions of the Einstein 



equations with a positive cosmological constant in an 'initial value problem' spirit, cf. [I8j (see also 




In this paper, we are interested in the possible genericity of regular solutions defined by formal 
series expansions in the context of vacuum models in higher order gravity derived from the analytic 
lagrangian R + eR^ theory. In particular, we prove that the higher order gravity field equations 
in vacuum admit a unique solution in the form of a regular formal power series expansion which 
contains the right number of free functions to qualify as a general solution of the system. This 
requires a careful function counting technique, and for this purpose it is necessary to develop a 
formulation of the theory as a system of evolution equations with constraints. 
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The plan of this paper is as follows. In the next Section, we write the higher order gravity 
field equations as a system of first order (ADM- type) evolution equations and constraints, as in 
an initial value formulation of the theory, and show that this system has the Cauchy-Kovalevski 
property. This allows us to count the true degrees of freedom necessary for any analytic solution to 
be a general one. In Section 3, we introduce a regular formal series representation of the metric and 
prove our main result about the existence and uniqueness of a regular generic perturbation in higher 
order gravity starting from given formal asymptotic data. Finally in the last Section, we present 
our conclusions and further discussion of these results. In Appendix A, we give more details about 
certain crucial steps in the proof of the Cauchy-Kovalevski property of our main system, while in 
Appendix B we present the full expressions of certain curvature components in terms of asymptotic 
data. 



2 Function counting 

In this Section, we are interested in counting the true degrees of freedom of the R + eR^ theory in 
vacuum. This requires a splitting formulation of the theory into evolution equations and constraints, 
which in turn relies on certain technical details that the resulting dynamical system has to satisfy 
for the whole scheme to be consistent. The main technical result is given at the end of this Section, 
Theorem 2.1. We then count the degrees of freedom of the theory. In the Appendix A, we give 
more details of the proof of the main theorem of this Section. 

We consider a spacetime (V,^) where V = M x 7W, with Ai being an orientable 3-manifold, 
the submanifolds Mt = {^i,^ M, t £ M, are spacelike and 5 is a Lorentzian metric, analytic and 
with signature (+, — , — , — jj. We take a Cauchy adapted frame Cj = (eo, e^) with tangent to the 
space slice A4t and eo orthogonal to it. The dual coframe 0* = {6^ = dt,6°' = dx°^ + (3°'dt), where 
the tangent vector /3" is the usual shift, leads to the standard general form of the metric 

ds^ = N^dt^ - gai3{t) (dx'^ + (3'^dt) (dxf^ + ^^dt^ . (2.1) 

Here A'^ is a positive function, the lapse, and we assume that all metrics gap{t) are complete 
Riemannian metrics. We could continue using a bar to denote spatial tensorj^, in a Cauchy adapted 
frame we would have gap = —gap and g°'^ = — g"^. However, we more frequently write = —dap- 
Below we shall deal exclusively with the simplest gauge choice, N = 1, (3 = 0, which means that 
goo = l,goa = 0, and local coordinates are adapted to the product structure on Aiti {x^ = {tjx")) 
(in this case, one sometimes speaks of a synchronous system of local coordinates). 



^Our conventions are those of [3]. 

^We generally follow closely the terminology of the fundamental treatise [1]. 
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Our starting point is the vacuum field equations of the higher order gravity theory derived from 
an analytic lagrangian f{R), that is the equations 

Lij = f'{R)Rij - ^f{R)gi, - V^VJf'{R) + giPgf'iR) = 0. (2.2) 

We shall be concerned below with the quadratic theory f{R) = R + eR^, for which, the field 
equations (|2.2p in a Cauchy adapted frame split as follows: 

Loo = (1 + 2eR)Roo - ^(1 + eR)R + 2eg''^VaVpR = 0, (2.3) 

Loa = (1 + 2ei?)i?oa - 2eVoV«i2 = 0, (2.4) 

L^p = (1 + 2eR)R^p - + ^R)R9ap - 2eV„V^i? + 2eg^pagR = 0. (2.5) 



The components of the Ricci tensor are 



Ri 



ldtK-\K^Kl (2.6) 



■00 - ~2 * ~ 4^a-^/3 
Roa = l{V0K^-VaK), (2.7) 

and also 

Rap = Pap + \dtK^p + \kK^p - \kIKp^, (2.8) 
where K = iiKap- Here, the extrinsic curvature is defined by the first variational equation 

dtlap = K^p, (2.9) 

and Pap denotes the three-dimensional Ricci tensor associated with ^ap- Further, we define the 
acceleration tensor D^p through the second variational equation 

dtKap = Dap, (2.10) 
and we also introduce the jerk tensor (3rd order derivatives) Wap through the jerk equation: 

dtDap = Wap. (2.11) 

The space tensors D^p and Wap are obviously symmetric and because of the fourth order nature 
of the higher order field equations, they play an important role in what follows. 

In terms of the acceleration tensor D^p, the Ricci tensor splittings given above become 

R,o = -^D + ^K^K^, (2.12) 
Roa = ^{VpK^-VaK), (2.13) 



Ra^ = Pap + \d^p + ^KK^p - ^K^K^^, (2.14) 
while for the scalar curvature we obtain, 

R = .p.D-^K' + ^K^K^,, (2.15) 

where P = ti Pa/3, D = tr D^is- If we substitute these forms into the identities (|2.3p . (|2.4p . we find 
the following equations, which, in obvious analogy with the situation in general relativity, we call 
constraints: 

Hamiltonian Constraint 

1„ l^.o 1 

8^ 



2^^NaVp{-P - ^K^ + ^K^'K^s - D)] = 0, (2.16) 
Momentum Constraint 

e[i-P - + \k'^K2 - D){VpKPa - V^K) - 

Va{-2dtP + KK^^K^s - 3K^K]Ki -KD + BK^^D^s - 2W)] = (2.17) 



The constraints show that the initial data (7^/3, i^a^, -Dq./?, PKva) cannot be chosen arbitrarily and 
must satisfy the equations ()2.16p and ()2.17p on each slice M.£.. 

Further, any initial data set {Ait, 1013, Da/3, W) must satisfy, apart from the equations (|2.9p 
(velocity equation), (|2.10p (acceleration equation), (j2.1ip (jerk equation), the following evolution 



■^It is to be noted that because of the structure of the last term in equation (|2.5p . the time derivative of the jerk 
tensor Wap appears only for its trace, whereas the time derivatives of the other components of Wajs do not appear. 
This is reminiscent of a similar fact in the FRW equations for universe with scale factor a and filled with a perfect 
fluid in general relativity. There, there appears no time derivative of the pressure in the Friedman equations. Of 
course, one does not count the pressure as an initial value in the evolution which is determined by, say, a,d,p as 
unknowns. Consequently, below we shall write {Mtjja^, Ka/3, Da/3 , W) for the initial data of the theory. 
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equation, called the snap equation: 

blK'^K^Dl + —DK'^^K^p - -D^ + IbD^^D^p - 3KW + 
l5K"^Waf,-6dt{dtP)- 

^r^VaVpi-P -D + Ik-^K^^ - \k^)] (2.18) 

The four evolution equations (|2.9|) . (|2.1U|) . (|2.1ip and (|2.18|) are the higher order gravity analogues 
of the ADM equations of general relativity, and together with the constraints (I2.16P and (I2.17P 
describe the time development (V, g) of any initial data set {M.t,^ai3, Kajj, Dap, W) in higher order 
gravity theories. To be able to complete the function counting argument of this Section, we need to 
prove that these equations are well defined in the sense that there is a well defined Cauchy problem, 
at least for the analytic case. 

To be concrete, one needs to prove that these equations are of the Cauchy-Kovalevski type, 
that is there are no time derivatives in the constraints and the derivatives of the unknowns are 
(through the evolution equations) analytic functions of the coordinates, the unknowns, and their 
first and second space derivatives. This is shown in the Appendix. From this result and the 
Cauchy-Kovalevski theorem (cf. [Ij, Appendix V), we are immediately led to the following result, 
local Cauchy problem (analytic case): 

Theorem 2.1 For N = 1, (3 = if we prescribe analytic initial data ("fais, Ka/s, Dai3,W) on some 
initial slice Mq; then there exists a neighborhood of Mq inMx M such that the evolution equations 
/i2.9\) . (2.10]) . \2.11\) and Ii2.18\) have an analytic solution in this neighborhood consistent with these 
data. This analytic solution is the development of the prescribed initial data on A4o if and only if 
these initial data satisfy the constraints. 

The last part of this theorem follows most easily from the conformal equivalence theorem of higher 
order gravity theories [20], that is working in the Einstein frame representation, and a theorem 
on symmetric hyperbolic systems, cf. [IJ, pp. 150-1, and also [21j. Indeed, in the Einstein frame 
representation, the theory is general relativity plus a self-interacting scalar field, and the system 
becomes one of the form 

i?o = -^dtK - ^K^K^ = 8vrfc(T° - ^T), (2.19) 
6 



Rl = \{VpKi - VaK) = SirkT^, (2.20) 
R^ = -P^ - ^dtiVlK^) = 87rk{T^ - ^jfT), (2.21) 

and the wave equation V*Vj</) — V'{(j)) = for (p having the particular scalar field potential given 
in [20]. This system is of the form given in Thm. 4.1 of [1], p. 150, from which the result follows. 

We are now ready to count the true degrees of freedom of the higher order gravity theories we 
consider. There are 19 relations, that is 18 by the three evolution equations ()2.9p . (|2.10p . (12. lip 
and one from (I2.18p . But we have the freedom to perform 4 diffeomorphism changes and we also 
have the 4 constraints ()2.16p and ()2.17p . Hence, the vacuum theory has 19 — 4 — 4 = 11 degrees 
of freedom. This in turn implies that any solution with 11 free functions has the same degree of 
generality with a general solution of the theory. 



3 Genericity of regularity 

In this Section we perform a perturbative analysis of our basic system of higher order equations, 
namely, the evolution equations (|2.9p . ()2.10p . (j2.1ip and (j2.18p together with the constraints (j2.16p 
and ()2.17p . In particular, we assume a regular formal series representation of the spatial metric of 
the form 

^ ^ _ ^(0) , (1) ^ , (2) ^2 , (3) 3 , (4) ,4 , . . . 

7a/3 — lap + 'aP ^ + 'a/3 + 'a/3 ^ + 'a/3 + ^'^■^1 

where the 7^'^^ , 7^^^^ , 7^^^ , 7^^^ , 7^"^^ , • • • are functions of the space coordinates. Because of the order 
of the higher order gravity equations, we shall be interested only in the part of the formal series 
shown, that is up to order four, and because of that we shall often drop the dots at the end of the 
various expressions to simplify the overall appearanccl. Thus before substitution to the evolution 
and constraint higher order equations, the expression ()3.ip contains 30 degrees of freedom. Note 
that setting 7^*^^ = dap and 7^'^-* = 0, n > 0, we have Minkowski space included here as an exact 
solution of the equations, and so our perturbation analysis covers also that case. 

The problem we are faced with in this Section is what the initial number of thirty free functions 
becomes after the imposition of the higher order evolution and constraint equations, that is how it 
finally compares with the 11 degrees of freedom that any general solution must possess as shown 
in the previous Section. Put it more precisely, given data aap , bap , Cap , dap , , analytic functions 



^Differentiation of such formal series witii respect to either space or the time variables is defined term by term, 
whereas multiplication of two such expressions results when the various terms are multiplied and terms of same 
powers of t are taken together. 
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of the space coordinates, such that the coefficients 'y'^^ , fJ- = 0, • • • 4, are prescribed, 

(0) (1) u (2) (3) , (4) „x 



how many of these data are truly independent when (j3.ip is taken to be a possible solution of the 
evolution equations ([23]), (f2ll)]l . (|2lT]) and (I2l8l) together with the constraints (|2J6]l and (IZTTll ? 

To proceed, we shall need the formal expansion of the reciprocal tensor 7°'^, where ja/sl^'^ = 'ho- 
using this, the various coefficients 7^^-* /i = 0, • • • , 4, of t in the expansion 7"'^ = X^J^g 7^"^ "'^^^ 
are found to be: 

+ (^-e"^ + b'^^d^ - b'^^b^^ ^ + + ^7 - 4^5 + ^°^4^5^f - ^^5) (3-3) 

Note that a^pa^'^ = 52 and the indices of 60/3,0^/3,^0^,6^/3 are raised by a°^^ . For any tensor X, 
using the formal expansion (]3.ip . we can recursively calculate the coefficients in the expansion 

= 4? + ^ + 4? + 4? + 4? + • • • , (3.4) 

in particular we can write down a general iterated formula for the n-th order term, . For 
instance, for the extrinsic curvature writing 

and using the first variational equation (|2.9p , the coefficients i^^'^'' are given by the general recursive 
formal expression 

7S=7S<^V (3-6) 
In terms of the data a, 6, c, , d, e, we have explicitly, 

Kal3 = dtlali = bap + 2Co/3t + 3(io/3t^ + 46^/3*^, (3.7) 

and for the mixed components we obtain. 



= j'^^K^p = 6^ + {2c$ - b''n^p)t + {M% - 26"^c^/3 + b'^'b'lb^p - c"^6^/3)t' 
+ {Ae^l - 36"Td^/3 + 26"'^6]'c^/3 - 2c"^c^/3 - ^"^6^/3 + 6"^cj6^/3 - b'^^blbjb^p 
+ c-'b]b^p)t\ (3. 



Further, setting 7 = |7a/3| = — ff, the mean curvature, 

K = K^ = r^daap = dM^), (3.9) 



is given by the form 

K = b+{2c-bp^)t + {3d-3b'^c^i + b''^b}b^a)t^ 

+ (4e - 46^c?^ + 46"^c^^a - b^^b^bt^bsa - '^c'lci)t\ (3.10) 

For completeness, we also give the expressions of the coefHcients of the fully contravariant 

symbols, 

^ _ 2(6"T5^ _ _ + b^'^C^ - b^^^b^^g + c'''^b^)f 

+ b"^bi^bHP^-d'^bt^bl)t^. (3.11) 

Using these forms, we can find the various components of the acceleration and jerk tensors to the 
required order. We have that the perturbation of the acceleration tensor in terms of the prescribed 
data is given by the form: 

Da/3 = dtKo^p = 2cap + Gdapt + 12ea^t'^ + ■■■ . (3.12) 

Further we find, 

= ^''^D^p = 2c^ + 2(3d^ - b^c})t + 2(6e^ - 36^d^ + bp^^c} - c'^c})^, (3.13) 

and 

D = 2c + 2{M - b1ci)t + 2(6e - 36^d^ + 6^6^c- - c^cf )^^ (3.14) 
where the trace is given by 

D = D2 = T^dtK^p. 
For the fully contravariant components, we find, 

^ 2c'^^ + 2{M'^^ -b'^'^c<^^-b^^c1)t 

+ 2(66°^^ - 36"Td^ - 3b^^d1 - 2c°'^c^ + b'^^bjc^ + b^f^b^c] + b^'^b^c2)t'^. (3.15) 
Lastly, the jerk perturbation series is found to be 

Wa(i = dtD^p = 6daf} + 24ea^t + • • • , (3.16) 

so that 

W = 6d + 6{4e - b'Sd^)t + ■ ■ ■ (3.17) 
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where 



w = w^ = r^dtD^f 



The components of the Ricci curvature are more comphcated when expressed in terms of the 
asymptotic data a, b, c, d, e and we give them in the Appendix B. In terms of those and the various 
expressions obtained above, the hamiltonian constraint (j2.16p becomes 



Co = lp-h% + h' 



+ e 



32 



16 



32 



+ 2a^"V^V,P - 2a^-V/3V„(-2c + -hi^hj - -/) 



+ 



+ 2R^^\-c + hjl^Xp) - R^^'^R^^'^ - 2a"'^V«V/3i?(^) + 26;^VaV^i?(°)] 1 1 + • • • = 0.(3.18) 



From (j2.17p . we calculate the momentum constraint in the form 
1 



Cn 



+ 6 [2fl(0) ((V^c^ - V,c) - ^V^ib^^bl) + lv.{b%)^ + R^'Hv^b^^ - V^b) 
- 4V„ii(2) + 5/3v^^(i) + (2c^ _ 6^62)V/3i?(°)] } t + • • • = 0. (3.19) 
Finally, the snap equation (|2.18p gives: 
lp + 2cj6^5^ + l.^ 

+ e |2[-P + (-2c + ^b'^b] - h^)] [-P + (-C + ift^^^ - h^)] -l[-P+ (-2c + ^b'^b] - \b')]' 

+ 6[-d^P - 3bd - 2^ - 24e + 216^4 + lOcJc^ - mibi^cj + ^6^625^6^ + 36]c^6 

+ 26]6^c - b1b%b - \{b]bi^f] + 4a"^[V^V„P - V^V„(-2c + -^bib% - \b^)]^ + • • • = (3.20) 

We realize that because of the presence of the e parameter in the theory, the resulting perturbation 
series become considerably more complicated in higher orders, due to terms of the order 0{et). To 
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simplify our further work, we find it convenient to use the fohowing notation: 

Co = A[J + (A Jge + (At)°t + (Arf)l]et + • • • 
C„ = A° +(Ae)°e+(A^)°^ + (Arf)°e^ + ••• 
L^ = A^ + (A,)^e + --- 
snap terms = A + (Ae)e + --- 

where the indices of the A-coefficients indicate the powers of t, e that appear in the series. For 
instance, the 0{t) term in the momentum constraint ()3.19p is denoted by 

We are now ready to decide how much the imposition of the field equations (j2.9p . (j2.10p . (j2.1ip 
and (|2.18p together with the constraints (j2.16p and (j2.17p . restricts the number of free functions 
in the data (j3.2p in the perturbation series (j3.ip . Using the field equations and calculating the 
various relations that appear at each order, we are led to the following relations between the data 
a, 6, c, d, e: From Aq, A, we get one relation, namely, 

c = \b^b^p, (3.21) 

from A^ we obtain three more relations, 

= V„6, (3.22) 

whereas from A, Aq and Eq. (|3.2ip we get the following six relations, 

-P^-^bP^b+h^^bl-ci = 0. (3.23) 

From A and Eq. (fOT]) . the (A,)^ term gives one more relation 

+ 6bd + 48e - m^d^ - 20c^c'^ + 38b^b^c} - 9b^blb2bf^ - Gb^c'^b + 2b]^b^b'^b = 0, (3.24) 

and from Eq. (|3.2ip . the {Xt)o term gives the relation 

26c + 2b^^c'^ - bib'^bl = 0. (3.25) 

From Eq. (I3.2ip and A, the term gives three more relations 

V„(6c - 6d + 56^c^ - hy^bt^bl) = 0. (3.26) 

Further, we note that the term (A^t)!] is identically zero due to the term A and the equations 
(|3.2ip . (|3.26p . The term (Af)[^ is identically zero due to the term A, the equations ()3.2ip . (j3.22p . 
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mjTI^ and the fact that V^P^ = iV^P. The term (Aet)° is identically zero due to the equations 
(|3.'2ip . (I3.22|) and the term {Xt)a- "^^^ terms {\e)o and are identically zero due to the equations 
([320, (13231) • 

Hence, in total we find that the imposition of the field equations leads to 15 relations between 
the 30 functions of the perturbation metric (j3.ip . that is we are left with 15 free functions. Taking 
into account the freedom we have in performing 4 diffeomorphism changes, we finally conclude that 
there are in total 11 free functions in the solution (|3.ip . This means that the regular solution (|3.ip 
corresponds to a general solution of the problem. Put it differently, 'regularity is a generic feature 
of the R + eR^ theory in vacuum'. 

We now ask: Out of the 30 different functions a, b, c, d, e, which eleven of those are we to choose 
to use as our initial data? We have shown in this Section that the vacuum higher order gravity 
equations 1^3^, ([2T0I1 . ([CT]) and (|2T8]) together with the constraints ([236]) and dST]), admit a 
regular formal series expansion of the form ()3.ip as a general solution requiring 11 smooth initial 
data. If we prescribe the thirty data 

aal3, ba^, Cap, daj3, e^^, (3.27) 

initially, we still have the freedom to fix 19 of them. We choose to leave the six components of 
the metric Ua/s free, and we choose the four symmetric space tensors ba/s^Caf^^daf^ and to be 
traceless with respect to Oq^. Then we proceed to count the number of free functions in several 
steps. In the first step, using the trace condition on 6q,^, 6 = 0, ()3.22p gives 

= 0. (3.28) 

Hence only two of the six components of ba/s can be free. Secondly, if c^^ has zero trace, then 
Eq. (|3.2ip fixes one more component of bap, thus we end up with 7 free functions (6 from a^^ and 
one more from 5^^). Then, as step three, we use the 6 relations in ()3.23p to completely fix the 
remaining 6 components of bap and c^^, namely, they fix the one remaining component of bap and 
the 5 ones of c^^ from step 2. Step four consists in taking into account the five remaining relations 
(I3.24p . (l3.25p and (I3.26P : We find that only seven from the twelve remaining components of the 
tensors da/s and e^/? can be independent. But since we took da/3 and Ca/s to be traceless, as a last 
step in our function counting, we find that together they can have 5 free, independent components 
at most. Summing up the free functions we found up to now, we end up with 

^ + ^ + ^ + ^ = 11 (3.29) 

from a^^ from b^^ from c^^ from d^^is and e^^ 
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suitable free data as required for the solution to be a general one. We thus arrive at the following 
result which summarizes what we have shown in this Section, and generalizes a theorem of Rendall 
|18j for higher order gravity theories that derive from the lagrangian R + eE? . 

Theorem 3.1 Let a^^ he a smooth Riemannian metric , bapiCap^dap and e^^ he symmetric smooth 
tensor fields which are traceless with respect to the metric a^/s, i-e., they satisfy b = c = d = e = 0. 
Then there exists a formal power series expansion solution of the vacuum higher order gravity 
equations of the form (3. 1\) such that: 

1. R is unique 

2. The coefficients ^('^)'^^ are all smooth 

3. R holds that 7^°j = a^p and 7^^j = 6^/3, 1^1 = ^ap, 1^} = dap and = e^p. 

In the course of the proof of this result, uniqueness followed because all coefficients were found 
recursively, while smoothness follows because in no step of the proof did we found it necessary to 
lower the C°° assumption. We also note that 6^/3 and Cap are necessarily transverse with respect to 

4 Discussion 

In this paper, we have treated the problem of the existence of generic perturbations of the regular 
state in higher order gravity in vacuum that derives from the lagrangian R + eR^ . We have shown 
that there is a regular state of the theory in the form of a formal series expansion having the same 
number of free functions as those required for a general solution of the theory. This means that 
there exists an open set in the space of initial data of the theory that leads to a regular solution 
having the correct number of free functions to qualify as a general solution. 

To achieve this, we have shown that there exists a first order formulation of the theory with the 
Cauchy-Kovalevski property. This formulation of the quadratic theory R + eR?' evolves an initial 
data set {M. , , K^p , D^p , W) through a set of the four evolution equations (j2.9p , (|2.10[) , (j2.1ip 
and (I2.18P and the two constraints ()2.16p and ()2.17p . and builds the time development (y,g)- What 
we have proved is that if we start with an initial data set in which the metric has the asymptotic 
form ()3.ip and evolve, then we can build an asymptotic development in the form of a formal series 
expansion which satisfies the evolution and constraint equations and has the same number of free 
functions as those of a general solution of the theory. In other words, we have shown that regularity 
is a generic feature of the R + eR^ theory under the assumption of analyticity. 
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The results of this paper provide the necessary background for various further investigations 
that we carry out currently. For example, elsewhere we plan to examine the generic perturbation 
problem of the well known radiation solution of the quadratic theory and extend the results of 
this paper to any perfect fluid spacetime with equation of state p = wp in higher order gravity. It 
is also interesting to further compare these radiation perturbations with the situation in vacuum 
in the context of higher order gravity. We also wish to extend our present results to the case of 
arbitrary lapse and shift, and to consider the problem of the present paper in the conformal frame. 
We know that in the case of a positive cosmological constant and a general perfect fluid source, 
similar results to those of this paper hold, cf. |18j . However, the scalar field case, especially with 
the potential of the Einstein frame representation of an f{R) theory, is to our knowledge an open 
problem. 
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Appendix A: Proof of the Cauchy-Kovalevski property 

In this Appendix we present details of the proof that the four evolution equations ()2.9p . ()2.10p . 
(j2.1ip and (|2.18p constitute a Cauchy-Kovalevski type system. Indeed, the only 'dangerous' terms 
present in these equations are the three terms 



P, dtidtP), V^Vp{-P-D + -^K^^K. 
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(4.1) 



present in the constraints and in the snap equation (j2.18p . 

The spatial connection coefficients are given by the obvious formula 



(4.2) 



We then have 





(4.3) 



and 



^t^ aB — 



i(L»^^ - 2i^^''/^^)(9^7,, + a„7;3, - a,7«/3 



) 



+ 



(4.4) 
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We also set 

= ^{dp^av + dalPr, - dr,^ap), (4.5) 

SO that = and further we set 

Zr)aP = -{dpKar, + daKp^ - drfK^^), (4.6) 

and 

Hr)ap = -^{d^Dan + daD^^ - djjDa/s)- (4.7) 
It then follows that the time derivatives of these 'fully covariant symbols' satisfy 







(4.8) 






(4.9) 




= H^a^. 


(4.10) 



Hence we conclude that the first time derivatives of the spatial connection coefficients depend 
only of (7a/3, -f^a/?) and their first spatial derivatives, and the second time derivatives of the spatial 
connection coefficients depend only of {"fapj K^js, D^js) and their first spatial derivatives. 
Now, the spatial Ricci tensor is given by 

PaP = - 9p^'a, + Kp^U - Ke^^, (4-11) 

and so its time derivative is calculated to be 

dtPap = d^^idtKp) - Mdtr^^) + dtKp^U + KpdtTU - dtF^^.F'^^^ - F^^A^^.- (4.12) 
Then, for the spatial scalar curvature, 

P = 7''^Pap, (4.13) 

we find that 

dtP = -K'^^Pap + r^dtPap, (4.14) 
and therefore the first of the dangerous terms finally reads: 

dtP = K<^^i-di,r'F,^p - r'di^F,^^ + dp^'^'F.ai, + I'^'dpF.a^^ - T^'T^T.^^r^^^ + 7'^''7<r^«.r^^^) 
+ 7"^[-5^K'^T,„;3 - K^^'di.F.^p + 5^7^'^6a/3 + I'^d^Z,^^ - dpKi^F,^^ - K^^'d/^F.^f, 

+ 9^7^^Z,,^ + 7^^a^Z,„^ 

+ 7^''r^/.c(-^'''r,„^ + i^'z^ap) + r^F^api-K^^'F^^e + Y'^Zr,^,) 

- l'^F^P^.i-KI'^Fr,ae+r''Zr,ae)-r%aei-K''^F^P^ + j'^Z^^^)]. (4.15) 
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This means that the first dangerous term is 'purely spatial'. This result also implies that the third 
dangerous term is also purely spatial, for it is calculated to be of the form 

- ^-K'^^K^K^p-^-KD+^-KK'-^K^p), (4.16) 

that is trouble could only had arisen from the first dangerous term, which however as we showed 
above is purely spatial. 
Lastly, since 

dt{dtP) = -{D^^ - 2K^^K^)P^p - 2K-^dtPap + 7"^^t(9tPa/3), (4.17) 

and 

d^p^p = 9^(9|r^^)-5^(5|r(i^) + a|r^^r;., + r^^5|r^,-a|r(i,r^^ 

- afr^^r(i, + 25tr^^5tr^,-2atr(i,atr^^, (4.i8) 

we find that the second dangerous term depends on (7^^, Kajj,Dajj) and its first and second spatial 
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derivatives, namely, it has the form: 

- ^<V^fi^{-K^^V^^, + 7^'' - 7^^r^,,(-K^''r,^^ + 7^''^^/3m)] 

+ 2dfiK^'Z,^^ + 2i^^^9/3Z,„^ - dfi^^^H,^^ - r'dpH,^^, 
+ 7C'?[-(I)A«^ - 2K>'^^)T,^p - 2K>''Z,o.p + r'H,o.p]^r,K 



+ _ 2K'^Kl])Tr,^, - 2K'^Z^^, + f'i/^^elrca 



/3 



+ 2{-K>''r,ap + 7'''^.a^)(-^^''r^MC + T^^'^wc) 

- 2(-i^^''r,,, + 7'^''Z^„,)(-i^<r^^^ + 7<Z^;3^)}. (4.19) 
This completes the proof that the evolution equations are a Cauchy-Kovalevski system. 

Apppendix B: Ricci curvarure in terms of the data a, 6, c, e 

We give here the various components of the Ricci curvature and the space-space components of the 
field equation ()2.5p in terms of the data o, b, c, d, e. We have: 

+ (-6e + - ^6^&^c^ + lb^,bZb2b$ + 2c^c^) ^ • • • , (4.20) 



17 



+ 

+ 

+ 
+ 



(4.21) 



+ [-Ui + 26?c^ - hlhi^bl + c^^bl - hie + \bt^blbi - \cih + \bP^blb^ t 
+ (-6ef + - 36?6^c2 + Sc^c^ + \d% - h^d^bl + ^6? 6f 6^62 - ^cf 6^6^ 

- \dib + ^6^c2:6 - -fi%b + i^fe^^ - c^c + \cibt,b] + i6^62:c - \bP^bl 



+ ••• , 



(4.22) 



where Pq^ the Ricci tensor associated with 7q,^. The scalar curvature becomes an expression of the 
form, 



+ 



(4.23) 



exphcitly we have, 



4 4 " 

= (-P - 2c + ^6^6^ - ^b'^ + (^-6c - 6d + 56^c^ - ^b^b^b} + 6^6^ i 
+ (^-^M -c'- 12e + + 5c^c^ - ^ 6^ 6^c^ + + ^^'^c^?' 

+ bib'^pc-\bXb^b-l{bib-0)Af + .... 



(4.24) 
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Using these forms, the space-space components of the field equations in a Cauchy adapted frame, 
()2.5p . become: 

+ e [n-P + (-2c + -fi] - -/)] [-PP^ + {-ci + \b^^hl - \hib)] 

- \[-P + (-2c + \ht^h] - \h')]Hi - 2a^^[V,V,P - (-2c + 6| - \h%^,,] 

+ 2[-^P - 36d - 2c2 - 24e + 216^4 + lOcJc^ - 195|6^c7 + ^ftlJ^J^^b^ + Sft^cl^ft 

+ 26]6?c - h]bihl^h - \{qh'^fn + 2a'<'[V,V,P - V,V,{-2c + - \h^)]5i^ 

+ • • • = 0. (4.25) 
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